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Abstract 

A consistent Hamiltonian formulation of the teleparallel equivalent 
of general relativity (TEGR) requires the theory to be invariant under 
the global SO (3) symmetry group, which acts on orthonormal triads 
in three-dimensional spacelike hypersurfaces. In the TEGR it is possi- 
ble to make definite statements about the energy of the gravitational 
field. In this geometrical framework two sets of triads related by a 
local SO(3) transformation yield different descriptions of the gravi- 
tational energy. Here we consider the problem of assigning a unique 
set of triads to the metric tensor restricted to the three-dimensional 
hypersurface. The analysis is carried out in the context of Bondi's 
radiating metric. A simple and original expression for Bondi's news 
function is obtained, which allows us to carry out numerical calcu- 
lations and verify that a triad with a specific asymptotic behaviour 
yields the minimum gravitational energy for a fixed space volume con- 
taining the radiating source. This result supports the conjecture that 
the requirement of a minimum gravitational energy for a given space 
volume singles out uniquely the correct set of orthonormal triads. 
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I. Introduction 

General covariance and the notion of absolute parallelism of vector fields in 
space-time are the two major geometrical concepts of teleparallel theories 
of gravity. M0ller|J was probably the first to make use of the concept of 
absolute parallelism to put forward an alternative geometrical framework 
for general relativity that could satisfactorily address the problem of the 
definition of the gravitational energy. 

The idea of absolute parallelism can be established by considering a space- 
time vector field and a set of orthonormal tetrad fields e a ^(x). At the 
space-time point x x the tetrad components of the vector field are given by 
V a (x) = e a ^(x)V^(x), and at x x + dx x by V a (x + dx) = e a ^(x + dx)V tM (x + 
dx) = V a (x) + DV a (x), where DV a (x) = e a ^(V \V^)dx x . The covariant 
derivative V is constructed out of the connection 

The vector field V^(x) is said to be autoparallel if its tetrad components at 
distant points coincide. Thus V^(x) is autoparallel if V ' vanishes. There- 
fore connection (1) defines a condition for absolute parallelism, or teleparal- 
lelism, in space-time. Such connection only makes sense if the tetrad field 
transforms under the global SO(3,l) group. 

A gravity theory based on (1) obviously depart from the Riemannian 
geometry because the curvature tensor constructed out of it vanishes identi- 
cally. Inspite of this fact, there does exist a theory based on (1) that describes 
the dynamics of the gravitational field in agreement with Einstein's general 
relativity. M0ller called such theory the "tetrad theory of gravity" , but for a 
long time it has been known as the teleparallel equivalent of general relativ- 
ity (TEGR) |fj, 0] . The Lagrangian density for this alternative description of 
general relativity is constructed by means of a quadratic combination of the 
torsion tensor T a ^ = d^e a v — d v e a M , which is related to the antisymmetric 
part of connection (1). 

It is possible to establish a theory for the gravitational field directly from 
(1). However, in order to make contact with more recent analysis, it is in- 
structive to approach the TEGR firstly considering it with a local SO(3,l) 
symmetry in the Lagrangian context. Eventually we will return to the geo- 
metrical framework determined by (1). 
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Although not extensively investigated in the literature, the TEGR in 
Lagrangian form has been considered as a viable formulation of the gravi- 
tational dynamics 0, [5], ^, M, [| inspite of troubles that may spoil the initial 
value problem|], pLOfl . Such problems take place if the Lagrangian density of 
the TEGR is formulated with a local SO(3,l) symmetry, in which case the 
torsion tensor is defined by T a ^ = d^e a v — d u e a ^+uj^ a b e b u — lo u a \ } e h M . The 
spin connection u^ ab is totally independent of e a M and satisfies the condition 
of vanishing curvature: R a b^u(^) = 0. 

In order to solve problems with respect to the initial value problem, the 
Hamiltonian formulation of the TEGR with the local SO (3,1) symmetry was 



considered [^TJ. By working out the constraint algebra one ultimately con- 
cludes that in order to formulate consistently the theory in terms of first class 
constraints it is mandatory to break the local SO (3) symmetry of the action 
in Hamiltonian form. Therefore a consistent Hamiltonian formulation of the 
theory displays invariance under the global SO (3) group that acts on triads 
restricted to three-dimensional spacelike hypersurfaces. 

The major feature of the Hamiltonian formulation of the TEGR is that 
the integral form of the Hamiltonian constraint equation can be written as 



an energy equation of the type C = H — E = 0[ 12 , 13, 14 1. The reason 
for this follows from the fact that the Hamiltonian constraint contains a 
scalar density in the form of a total divergence, whose integral over the 



whole three-dimensional hyper surf ace yields the ADM energy 1 15]. As we 



will argue, the integral of this scalar density over finite volumes of the three- 
dimensional hypersurface yields a natural definition for the gravitational field 
energy. Therefore in the TEGR one can make definite statements about the 
localizability of the gravitational energy, inspite of claims according to which 
the latter is not localizable. In fact the very concept of a black hole lends 
support to the idea that the gravitational energy is localizable, since there is 
no process by means of which the gravitational mass inside a black hole can 
be made to vanish. A remarkable application of the energy expression of the 
TEGR has been made in the evaluation of the irreducible mass of a rotating 
black holedTB . 



The TEGR exhibits two specific properties: the emergence of a possible 
definition for the gravitational field energy and the global SO (3) symmetry of 
the theory. We believe that these two features are intimately related. Since 
the symmetry of the theory is global, triads related by a local SO (3) trans- 
formation are inequivalent and a priori we have no means to select the one 
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that actually describes the spacelike hypersurface. We conjecture that the 
requirement of a minimum gravitational energy for a given space volume is 
one condition that singles out uniquely the correct set of triads. In this pa- 
per we investigate this conjecture in the realistic context of Bondi's radiating 



metric [I7J. This conjecture was already put forward in a previous investiga- 



tion of Bondi's energy in the framework of the TEGR[18|. An additional 



essential requirement for a consistent expression of the gravitational energy 



is the boundary conditions on the triads. It was noted |12| that the ADM en- 
ergy is obtained from the energy expression of the TEGR if the asymptotic 
behaviour of the triads at spacelike infinity is given by 

e (i)j ~ Vij + gM~) > ( 2 ) 

irrespective of any symmetry of the tensor hij. 

A second possible, independent condition for assigning a set of triads to 
a given three-dimensional metric tensor amounts to requiring a symmetric 
tensor = hji in the asymptotic expansion of e^j. We will prove that 
this symmetry condition uniquely associates e^y to the metric tensor of the 
spacelike hypersurface of Bondi's metric. The unique character of such triads 
strongly supports this second conjecture. 

The two conjectures above are not mutually excluding. In this paper 
we argue that the set of triads that yield the minimum gravitational en- 
ergy within a space volume containing the radiating source is the one whose 
asymptotic behaviour is determined by the symmetry condition = hji. 
We will show that this fact is indeed verified by analyzing several configura- 
tions for the triads. Unfortunately we have not found it possible to prove on 
general grounds that the "symmetrized" triad yields the minimum energy. 

In order to obtain numerical values for the gravitational energy for a 
finite three-dimensional volume of Bondi's space-time we need an explicit 
expression of the news function. However the existing expressions in the 
literature are not suitable for our purposes. In particular, the expression of 



the news function given by Hobill|T9[ (to be presented ahead) is very intricate 
to the extent of not allowing a computer evaluation of numerical values of 
the gravitational energy. Therefore we have obtained an original and simpler 
expression for the news function that: (i) satisfies all necessary regularity 
conditions related to the axial symmetry of the system, (ii) makes the initial 
and final states (space-times) described by Bondi's metric to be nonradiative 
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and (iii) yields an expression for the energy density that can be numerically 
integrated. 

In Section II we describe the TEGR in Lagrangian and Hamiltonian for- 
mulations, show the emergence of the definition of the gravitional energy 
and further discuss the troubles that arise in the initial value problem of the 
theory if it is constructed with a local SO (3,1) symmetry. In section III we 
present Bondi's radiating metric and three expressions for triads restricted to 
the three-dimensional spacelike hypersurface. In this section we also prove 
that the symmetry condition hij = hji uniquely associates a set of triads 
(whose asymptotic behaviour is given by (2)) with the metric tensor for the 
spacelike section of asymptotically flat space-times. The news function and 
the related mass aspect that will be needed for the calculations of the gravi- 
tational energy are obtained in section IV. In section V we carry out several 
calculations that lead to the main conclusion regarding the selection of triads. 

Notation: spacetime indices fx, v, ... and local Lorentz indices a,b, ... run from 
to 3. In the 3+1 decomposition latin indices from the middle of the alphabet 
indicate space indices according to = 0,z, o = (0), (i). The tetrad field 
e a n and the spin connection u^ ab yield the usual definitions of the torsion 
and curvature tensors: R a b/il/ = d il u v a b + uj^ a c uj u c b - T a ^ = d^e a v + 
Wjj^e^ — .... The flat space-time metric is fixed by ?7(o)(o) — — 1- 



II. The Lagrangian and Hamiltonian formulations of the TEGR 

The Lagrangian density of the TEGR in empty space-time, displaying a 
local SO (3,1) symmetry, is given by 

L(e,u,\) = -ke( l -T ahc T abc + l -T ahc T bac - T a T a ) + e\ ah ^ R ab ^(u) . (3) 

where k = j^q, G is the gravitational constant; e = det(e a M ), \ ab ^ v are 
Lagrange multipliers and T a is the trace of the torsion tensor defined by 
T a = T b ba . The tetrad field e a ^ and the spin connection u^ ab are completely 
independent field variables. The latter is enforced to satisfy the condition of 
zero curvature. Therefore this Lagrangian formulation is in no way similar 
to the usual Palatini formulation, in which the spin connection is related to 
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the tetrad field via field equations. Later on we will introduce the tensor Ti abc 
defined by 



rpabcrp | rpabcrp rparp rpabc^ 

1 abc ~r — J- J-bac ~ 1 1 a = 1 ^abc ■ 

The equivalence of the TEGR with Einstein's general relativity is based 
on the identity 



eR(e,u) = eR(e) + e{-T abc T abc + l -T abc T bac - T a T a ) - 2d,(eT») , 

which is obtained by just substituting the arbitrary spin connection uj^ = 
°^afe(e) + K^ a b in the scalar curvature tensor R(e,u) in the left hand side; 
°^Vfc(e) is the Levi-Civita connection and K^ ab = \e a x e b v {T x ^ + T uXfl - 
T^vx) is the contorsion tensor. The vanishing of R a bfJiU (uj), which is one of 
the field equations derived from (3), implies the equivalence of the scalar 
curvature R(e), constructed out of e a M only, and the quadratic combination 
of the torsion tensor. It also ensures that the field equation arising from the 
variation of L with respect to e a M is strictly equivalent to Einstein's equations 
in tetrad form. Let j-^- = denote the field equations satisfied by e afl . It 
can be shown by explicit calculations that 

<)L 1 e{R a ,{e) - l -e a ,R{e)} . (4) 



5e a » 2 L *** w 2 
We refer the reader to Ref. ]TTJ for additional details. 



Throughout this section we will be interested in asymptoticaly flat space- 
times. The Hamiltonian formulation of the TEGR can be successfully im- 
plemented if we fix the gauge u 0ab = from the outset, since in this case the 



constraints (to be shown below) constitute a first class set[|TTJ. The condition 
^oab — is achieved by breaking the local Lorentz symmetry of (3). We still 
make use of the residual time independent gauge symmetry to fix the usual 
time gauge condition e^) = emv = 0. Because of u 0ab = 0, H does not 
depend on p kab ^ the momentum canonically conjugated to uikab- Therefore 
arbitrary variations of L = pq — H with respect to P kab yields Ukab = 0. Thus 
in view of u>o a b = 0, ujkab drops out from our considerations. 

As a consequence of the above gauge fixing the canonical action integral 
obtained from (3) becomesjll 
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Atl = J d A x{U^ k e m - H} , 



(5) 



H — NC + N*C t + S mn n m " + J-d k (NeT k ) + d k {W k N ) . (6) 

07TG 

N and N l are the lapse and shift functions, and E mn = — S„ m are Lagrange 
multipliers. The constraints are defined by 

C = d 3 (2keT>) - keE^Tuj - 7^( n %* " \^) > ( 7 «) 

C k = -e {3)k dM 3)l - n U) %)ik , (7b) 
where e = det(e (j)k ), T l = g' lk e i3)l T (j) i k , T (i)jk = dje (i)k - d k e (i)j , and 

Y^ijk _ }_(rpijk _|_ rpjik _ rpkij^ _|_ }i^ikrpj _ 

We remark that (5) and (6) are now invariant under global SO (3) and 
general coordinate transformations. Therefore the torsion tensor restricted 
to the three-dimensional spacelike hypersurface is ultimately related to the 
antisymmetric part of the spatial components of (1). Had we dispensed with 
the connection from the outset we would have arrived at precisely the 
canonical formulation determined by (5), (6) and (7). Such connection has 
been considered in previous investigations of teleparallel theories, but in fact 
it is eventually unnecessary for the establishment of the theory. 

If we assume the asymptotic behaviour 

e 0> ~ Vjk + \h jk (^) (2) 
for r — > oo, then in view of the relation 

y* d 3 xdj(eT J ) = J s dS k (dih ik - d k hii) = E ADM (8a) 

where the surface integral is evaluated for r — > oo, the integral form of the 
Hamiltonian constraint C = may be rewritten as 
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I d 3 x{keZ kij T kij + ^(n%* - -IT 2 )} = E ADM . (86) 

The integration is over the whole three-dimensional space. Given that dj (eT- 7 ) 
is a scalar density, from (8a,b) we define the gravitational energy density en- 
closed by a volume V of the space as 



It must be noted that E depends only on the triads e(k)i restricted to a three- 
dimensional spacelike hypersurface; the inverse quantities can be written 
in terms of e^v. From the right hand side of equation (4) we observe that 
the dynamics of the triads does not depend on uj^. Therefore E given above 
does not depend on the fixation of any gauge for uj^ ab . The reference space 



which defines the zero of gravitational energy has been defined in ref. 
We briefly remark that the differences between (9) and Miller's expression 
for the gravitational energy have been thoroughly discussed in ref. 



We make now the important assumption that the general form of the 
canonical structure of the TEGR is the same for any class of space-times, 
irrespective of the peculiarities of the latter (for the de Sitter space|2D[, for 
example, there is an additional term in the Hamiltonian constraint C). There- 
fore we assert that the integral form of the Hamiltonian constraint equation 
can be written as C = H — E = for any space-time, and that (9) represents 
the gravitational energy for arbitrary space-times with any topology. 

We recall finally that Miiller-Hoissen and Nitsch[jTU[ and Kopczyhski|| 
have shown that in general the theory defined by (3) faces difficulties with 
respect to the Cauchy problem. They have shown that in general six compo- 
nents of the torsion tensor are not determined from the evolution of the initial 
data. On the other hand, the constraints of the theory constitute a first class 
set provided we fix the six quantities uo a b = before varying the action ^T^ . 
This condition is mandatory and does not merely represent one particular 
gauge fixing of the theory. Since the fixing of u>o a b yields a well defined the- 
ory with first class constraints, we cannot assert that the field configurations 
of the latter are gauge equivalent to configurations whose time evolution is 
not precisely determined. The requirement of local SO(3,l) symmetry plus 
the addition of X ab ^R abflu (u) in (3) has the ultimate effect of discarding the 
connection. 
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Constant rotations constitute a basic feature of the teleparallel geometry. 
According to M0ller0, in the framework of the abolute parallelism tetrad 
fields, together with the boundary conditions, uniquely determine a tetrad 
lattice, apart from an arbitrary constant rotation of the tetrads in the lattice. 



III. Bondi's radiating metric and the associated triads. 

Bondi's metric describes the asymptotic form of a radiating solution of 
Einstein's equations. It is not an exact solution; it holds only in the asymp- 
totic region. In terms of radiation coordinates (u,r,9,<fi), where u is the 
retarded time and r is the luminosity distance, Bondi's metric is written 

as[ 



ds i = -lT e 213 -U 2 r 2 e 2 Adu 2 -2e 2f} dudr -2U r 2 e^dudO 



r 



+r 2 e 2 ~< d6 2 + e" 27 sin 2 6 dcf) 2 ) . (10) 



This metric tensor displays axial symmetry and reflection invariance. By 
requiring u = constant, (10) describes null hypersurfaces. Each null radial 
(light) ray is labelled by particular values of u, 9 and </>. At spacelike infinity 
u takes the standard form u = t — r. The four quantities appearing in (10), 
V, U, (3 and 7 are functions of u, r and 6. A more general form of this metric 
has been given by Sachsf21|j, who showed that the most general metric tensor 



describing asymptotically flat gravitational waves depends on six functions 
of the coordinates. 

The functions in (10) have the following asymptotic behaviour: 



V 2M 1 

9 
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U = - i (I + 2c cotf) + I (2d + 3c | + 4c 2 cotfl 

where M = M(u, 9) and ti = d(u, 9) are the mass aspect and the dipole 
aspect, respectively. From the function c(u,9) we define the news function 

dcjufi) 
du 

The functions U, V, (3 and 7 must satisfy regularity conditions along the 
z axis (9 = 0, 7r). We must require 

£/" 7 

V, 0, 



' ' sm0' sin 2 9 
to be regular functions of cos9 for 6* = 0, tt. The regularity conditions will be 
necessary for the construction of the news function, in section IV. 

The application of (9) to Bondi's metric requires transforming it to spher- 
ical coordinates (t, r, 9, </)) for which t = constant defines a space-like hyper- 
surface. Therefore we carry out a coordinate transformation such that the 
new timelike coordinate is given by t = u + r. We arrive at 

rfs 2 = _(V e 2P_ U 2 r 2 e 2 7 ) dt 2_ 2Ur 2 e 2'y dtdd 



dr dt 



+2 



e 2P[ V __ X \_ V 2 r 2 e 2 1 



+ 



,2/3, 



V 



2--) + U 2 r 2 e 2 ^ 
r 



dr 2 + 2U r 2 e 2l dr d9 + r 2 (e 2 ^d9 2 + e" 27 sin 2 9 djA . 

(11) 

Therefore the metric restricted to a three-dimensional spacelike hypersurface 
is given by 



ds 2 



V 



e 2P( 2 - -)+U 2 r 2 e 27 
r 



dr 2 + 2Ur 2 e 2 " 1 dr d9 



+r 2( e 2l d Q2 + e -2 1 sm 2 Qi 



(12) 



We recall that Goldberg [^2] and Papapetrou[|23[ have already considered 
Bondi's metric in cartesian coordinates. 
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The crucial point of the present investigation is the determination, in the 
framework of the TEGR, of the correct set of triads that lead to (12). If 
the metric tensor has only diagonal components, such as the metric tensor 
restricted to the three-dimensional spacelike section of Kerr's space-time, 
then the simplest construction has proven to be the correct one[16|. However, 



for metric tensors that contain off-diagonal terms, the determination of the 
unique triad is by no means a trivial procedure. In fact there is an infinity 
of triads that satisfy the boundary conditions determined by (2) and lead 
to (12). In Ref. fl8| two sets of triads that comply with (2) are presented. 
They are given by 



'(fc)i 



where 



/ A sinO coscf) + B cosO coscf) rC cosO cose, 
A sinO sincf) + B cosO sincfi rC cosO sine, 
\ A cosO — B sinO —rC sinO 



—rD sinO sine, 
rD sin9 coscf) 




(13) 



A = eM2 



V 



B = rUe 1 , 



C = e\ 



(14a) 
(Ub) 
(14c) 



D 



and 



(14d) 



where 



/ A' sinO coscf> rB' cosO coscf) + rC sin8 cose, 
A' sinO sincf) rB' cosO sincf) + rC sinO sine, 
\ A' cosQ -rB' sinO + rC cosO 



—rD' sinO sinCf 
rD' sin9 coscf) 




(15) 



A' 



2 - 



V 



+ U 2 r 2 e 2 ^ 



(16a) 
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(166) 

v /■ 

C = —Ure 2 \ (16c) 
A' K ' 

D' = e~ 7 . (16d) 

It is easy to see that both (13) and (15) yield the metric tensor (12) through 
the relation e^e^k = 9jk- They are related by a local SO (3) transformation. 

We have presented (13) and (15) because they are the simplest construc- 
tions that satisfy two basic requirements: (i) the triads must have the asymp- 
totic behaviour given by (2); (ii) by making the physical parameters of the 
metric vanish we must have T(}.)ij = everywhere. In the present case if we 
make M = d = c = both (13) and (15) acquire the form 

f sinO coscf) rcosO coscf) —r sinO sincj)\ 
£(k)i — sinO sincp r cosO sincf) r sinO coscf) . (17) 
\ cosO — r sinQ J 

In cartesian coordinates the expression above can be reduced to the diagonal 
form erk)i(x, y, z) = 5^. The requirement (ii) above is essentialy equivalent 
to the establishment of reference space triads, as discussed in [[L6| . A proper 
definition of gravitational energy requires the notion of reference space triads, 
which in the present case are given by (17). Note that by a suitable choice 
of a local SO(3) rotation we can make the flat space triads (17) satisfy the 
requirement (i), but not (ii). 

We proceed now to obtain the triads whose aymptotic expansion is given 
by (2) with the symmetry condition hij = hji. The procedure is the following. 
We consider, for instance, triads (13) and transform it to cartesian coordiates. 
Then we perform a local, asymptotic transformation 

e( k )i(t, x, y, z) = A (j) {k) e (j)i (t, x, y, z) , (18) 
where satisfies 

A (J) (fc) w 5 U) (fc) + {k) , (19a) 
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and W(j)(k) ~ 0(~) for r ^ oo (r = y/x' 2 + y' 2 + z 2 ). Transformation (19) 
preserves the asymptotic behaviour of the triads. By substituting (19) in 
(18) we find 

e(k)i = e(fc), + LO( k )i , 

from what follows 

hki = h ki + 2u ki , (20) 

where hki is given by the asymptotic expansion of (13) in cartesian coordi- 
nates. By requiring 

hki = h ik , (21) 
and making use of (19b) we find that 

^ki = -^{hki — hik) ■ (22) 
Substituting now (22) in (20) we arrive at 

hki = -Ahn + hik) = h(ki) ■ (23) 
Thus we ultimately obtain the symmetrized triads in cartesian coordinates: 

e( k )i ~ Vki + ^hi ■ (24) 

It must be noted that we arrive at a symmetrized triad irrespective of 
the triads we consider initially. The only requirement is that the unrotated 
triad must satisfy the asymptotic behaviour (2). Had we considered (15) we 
would arrive at the same result. In particular, from (22) we observe that if 
the initial triad is already symmetrized, then no local, asymptotic rotation 
is necessary. 

By transforming (24) into spherical coordinates t, r, 9, (/) (in which case 
the triads are no longer symmetrized), we finally arrive at 
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M f 

etui ~ (1 H )sin9cos<p cos9 coscj) , 

v > r 2r 

M f 

6(2)1 ~ (1 H )sin6 sintfi cos9 sinS , 

w r 2r 

M f 

e (3)i ~ (1 H )cos6' H sm^ , 

r 2r 

c f i 

6(1)2 ~ ^(1 H — )cos9 cos(p — —sin9 coscj) , 
r 2 

c / 
6(2)2 ~ r (l H — )cos6' sin0 — —sin6 sincf) , 

c f 

6(3)2 ~ —r(l H — )sin9 cos9 , 

r 2 



£(1)3 ~ — r(l )sin9 sincf) , 

6(2)3 ~ r (l )sin9cos(p , 



e (3 )3 w , (25) 



where / is given by 



Oc 

f = — + 2ccotg9 . 

We observe that by making M = c = triads (25) reduce to (17). We also 
note that (25), as well as (13) and (15), only make sense in the asymptotic 
region where Bondi's metric is valid. 

In view of (23) and (24) it is now easy to prove the uniqueness of the 
symmetrized triad. Suppose that the asymptotic behaviour of the metric 
tensor is given by 

9ij ™ Va + h*ij(-) . (26) 
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On the other hand by making use of (2) it follows from the relation = 
i e (k)j and from (23) that 

9ij « Vij + 7^( h ij + hji) = rjij + h(ij) . (27) 

Since h*j is unique, by comparing (26) and (27) we are led to conclude that 
there exists a unique symmetrized triad associated to the spatial section of 
an asymptotically flat metric tensor. 



IV. Construction of the news function and of the mass aspect 

In order to establish explicit expressions for the functions c(u, 0) and 
M(u,9), the suplementary field equations i?oo = R02 = for Bondi's metric 
are considered. In simplified form they are given by[|Hj| 



dM 



du 



dc\ 2 Id 



dd 



du) 
dM 



+ 



2du 



d 2 c dc 

w+3 -co tg e-2c 



du d9 



3 _|_ 4 ®° l Q _)_ ®° ®° 

dO du du dO du 



(28) 



(29) 



From (28) we observe that if J| 



0, the mass aspect M does not depend on 



u. 



For a family of null hypersurfaces Bondi's mass is defined by 



m\u) 



M(u,0) sinOdO . 



(30) 



It represents the mass of the system at the retarded time u. Multiplying 
both sides of (28) by sinO, integrating in and making use of the regularity 
conditions on the z axis stated in section III, we arrive at 



dm 
du 



dc\ 2 

— I sinO dO 
du) 



(31) 



which expresses the loss of mass. This equation was first obtained by Bondi 
et. al.[I7 . Therefore if the news function is nonvanishing, the mass of the 
system decreases in time. We remark that not only c(u, 0), but also the news 
function must be a regular function on the z axis. 
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The function c(u, 9) determines not only the loss of mass, but in fact it 
determines the whole structure of Bondi's metric, since via (28) and (29) 
it also determines the functions M, d, and all other functions that arise in 
the process of integration of the field equations. Therefore its construction 
deserves special attention. 



The news function proposed by Bondi et. al.[|T7]] is given by 

dc 00 

T- = fn{u)Pn(n) , 

du ^0 

where /i = cos9 . P n (fJ>) are the Legendre polynomials and f n (u) are func- 
tions to be determined. These functions become more and more intricate 
for increasing n, and eventually the above expression cannot be manipulated 
analytically. 

Bonnor and Rotenberg[^4j], Papapetrou|25| and Hallidy and Janis pjf 



have attempted at establishing an expression for c(u,9), but none of these 
proposals worked out satisfactorily, either because of the complexity of the 
structure of c(u, 9)^24^, or because the loss of mass does not occur for for 
a finite number of terms in the expansion in n[2"5|, Bonnor and Roten- 



berg's expression describes a radiative period between an initial nonradiative, 
static state and a final nonradiative, nonstatic state. The other approaches 
mentioned above describe a radiative period between static initial and final 
states. 

Hobill|]19[ has provided an expression for c(u, 9) that describes a radiative 



period between initial and final static states, and that leads to a loss of mass 
that is exactly equal to the total variation of the mass aspect during the 
period. Again the expression is not simple. It is given by 

( x _ K/(/i)e^ + m a (l - /i 2 )](l + W 2 )(l - /i 2 ) 
C[U,fl> / 2 (/i)e 2 ^ + (n-l)^ 2 -n^ + i ' W 

where n and r] are constants and m a is identified as one fourth of the total 
mass loss. It relates to the other constants according to 

18n 

r)m a , 



n + 1 
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/(/i) is an arbitrary function that must be everywhere regular and positive 
definite in the interval — 1 < /jl < 1. We note that c vanishes over the 
symmetry axis (/j, = ±1). The news function obtained from (32) reads 



dc 
du 



[e^fi^il + nfi 2 )- 1 + (1 
2m a (l-u 2 )f(u)e^ 



Z^ 2 )] 2 



-m b f 2 (u)e 2r,u 
1 + ny? 



m b (l - u 



(33) 



1 + nfi 2 

Expressions (32) and (33) lead, via integration of equation (28), to an ex- 
tremely complicated expression for the mass aspect: 

M(u,n) = M(-oo) -Am a 



1 + na 2 



+ (1 - ^ 



4m„(l 



2\3 



ll/U 



1=1 



(34) 



The intricate expressions for Hi(fj) are given in the Appendix of Ref. [1^ 



The problem with (32), (33) and (34) is that they are too complicated to 
yield numerical values for integrals (to be considered in the next section) 
that contain these expressions, even by means of computer calculations. 

Therefore we attempted at obtaining a simpler expression for c(u, fi). We 
damanded two conditions on c(w,/i). First, it must satisfy the regularity 
conditions on the z axis, which guarantees that the system is permanently 
isolated and leads to a well defined loss of mass. Second, that the initial and 
final states are nonradiative. However, as we will show, our expression leads 
to a nonstatic final state. It is given by 



ae nu {\ - fi 2 )F(fx) 



2nu 



+ 1 



(35) 



where n and a are constants (n^ 1 and a have dimension of length) and F([i) 
is a function that must be choosen such that c(l — /i 2 ) -1 = c(sin8)~ 2 is 
a regular function. Note that (35) vanishes for /i ± 1. The news function 
associated to (35) is given by 



dc 
du 



na{\ - y 2 )F{ii)e r - 



2nu 



+ 1 



1 - 



2e 

,2nu 



2/i u 



+ 1 



(36) 
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In the limit u — » ±00 we have 



dc 



lim c = lim — — = , 



(37) 



for any nonvanishing value of n. The property above is a necessary but not 
sufficient condition for having a static final state. Equation (37) does not 
determine whether the mass aspect M depends on 9, and so from (29) there 
is the possibility that d depends on u. The initial state is assumed to be of 
the Schwarzschild type. 

Integrating equation (28) from —00 to u, considering c(u, fx) given by 
(35), we obtain 



/u 
-ex 



2e 



3nu 



D 2nu 



+ 1 (e 2nu + If 



du 



+- 



1 ae r 



2 (e 2nu + 1 



1 - fx' fF" - 8/i(l - fi 2 )F' + A(3fi 2 - l)F(fx) 



+ A(/i) , (38) 



where F' — F" — and X(fx) is an integrating function that depends 
only on fx. 

The integral in (38) depends on the sign of n. For n<0we have 



M(u, fx) = M(-oo) + -na 2 {\ - fx 2 ) 2 F 2 (fx) 

6 



(e 2nu + If 



1 ae nu 



2 (e 2nu + 1) 
and for n > 0, 



1 - fx 2 ) 2 F" - 8fi(l - fx 2 )F' + 4(3fx 2 - l)F(fx) 



+ X(fx), (39) 



M(u, fx) = M(-oo) + na 2 (l - fi 2 ) 2 F 2 (fx) 



1 3e 4nu + 1 23 
6 (e 2nu + l) 3 + 48 



1 



ae 



2 (e 2nu + 1} 



1 - fx 2 ) 2 F" - 8/i(l - /x 2 )F' + 4(3^ 2 - l)F(/i) 



+ A(^) . (40) 
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By requiring the inital state to be static, expressions (39) and (40) must 
satisfy 



lim M = M(-oo) = M . 



(41) 



Applying the condition above to (39) we find that the function A(/x) vanishes, 
and therefore 



M(u, n) = M + -na 2 (l - /x 2 ) 2 F 2 (/i) 
o 



3e 



Anu 



(e 2nu + l) 3 



1 ae nu 



2 (e 2nu + 1) 



1 _ ^y F » _ 8/i (i _ ^) F > + 4(3^ _ i) j p( /U ) 



(42) 



for n < 0. Similarly, applying (41) to expression (40) we find that the 
integration function X(/j) is given by 



A( M ) = -(i + |)na 2 (l-MW) 
Substituting it back in (40) for n > we arrive at 



M(u, fi) = M + -na 2 (l - ^fF 2 ^) 



3e 4nu + 1 
(e 2nu + l) 3 



+ 



1 ae r - 



2 (e 2mt + 1) 



1 - /i 2 ) 2 F" - 8/x(l - fi 2 )F' + A(3fi 2 - l)F((i) 



(43) 



Let us now check the limiting value of (42) and (43) for u — > oo. Consid- 
ering first (42) we obtain 

lim M = M + -na 2 (l - ^ 2 ) 2 F 2 U) , 

where n < 0. It follows from the expression above that the total variation of 
the mass aspect AMt associated with (42) is given by 



AM T = ^na 2 (l- f i 2 ) 2 F 2 ( f i). 



(44) 
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The mass aspect M(u,/i) given by (43), for which n > 0, leads to a similar 
expression for A My: 



AM T = --na 2 {l - ^) 2 F 2 (^) . (45) 
6 



Hobill]T9j] obtained a result analogous to (44). However his expression, in 
the limit u — > oo, does not depend on fj,. Consequently the final state in his 
approach is static. In the present case, the mass aspect will still depend on 
/i in the limit u — > oo. In view of (29) this fact implies a final nonradiative, 
nonstatic state, which according to Hobill may describe a even more realistic 
situation. In any event expressions (35) and (36) describe an isolated physical 
system whose initial and final states are nonradiative, and whose loss of mass 
is precisely determined. 



V. Evaluation of the gravitational energy 

In this section we obtain numerical values for the gravitational energy 
given by expression (9). The relevance of such expression is that we can 
apply it to finite spacelike volumes. In the present case we will calculate 
the gravitational energy inside a large but finite surface of constant radius 
r , centered at the radiating source. Expression (9) will be evaluated as a 
surface integral in the asymptotic region where the components of Bondi's 
metric are precisely determined. It reduces to just one integral given by 

E = — [ dOdc^eT 1 , (46) 
8n Js 

where S is a surface of fixed radius r , assumed to be large as compared with 
the dimension of the source, and the determinant e is given by e = r 2 A sin9. 

We will consider triads (13), (15) and (25), for which the functions c(u, 9) 
and M(u, 6) are given by (35) and (42), respectively. Thus we take the mass 
aspect determined by the condition n < 0, only. Moreover we will consider 
four possibilities for the function F(n). Therefore we are effectively analyzing 
twelve distinct sets of triads. 



We will follow here the steps of section V of ||18|| . The energy expression 
that results from (13) and (15) have already been evaluated. They are given 
by0 
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E 1 = — F d6\ sinO 
4 Jo 



pi _|_ e -7 



i a 



(47) 



and 



sinO 



e^A' + e-^A'B' -2 + e~^^-C - BC' - Be~^ 

89 89 



-BcosO 



B' 



(48) 



together with definitions (14) and (16). E\ and E2 correspond to (13) and 
(15), respectively. Unfortunately the sign of the first term in the expansion 



of the function U(u,r,9) of the metric tensor (10), in Ref. []18l, is changed. 
Therefore expressions E\ and E 2 given in the latter reference must be cor- 
rected. Minor modifications (such as the modification of some numerical 
coefficients) are necessary. The correct expressions of E\ and Ei in terms of 
c(u, 9) and M(u, 9) are given by 



Ex = - r d9 sin9 M — F dO sin9 

2 Jo 4r Jo 



|) 2 + 4c(§)^H-4cW* 



1 r 8 
+— / d9 M — 
4r Jo 89 



fdc 

s i n ® tttt + 2c cot9 
\89 



(49) 



E 2 = - Fd9Msin9-— F d9 sin9 
2 Jo 4r Jo 



Ijr2 l/8c\ 2 A fdc\ n 



2\d9J 



2 2 „ f8M\fdc n 
+6 C 2 co^+^j(^- + 2 CC0 ^ 



1 r 71 " 

H / a!6>co,s6» 

4r ./o 



2c 



/dc 



4c 2 cot6> 



(50) 

It is clear that the total gravitational energy given by both (49) and (50), in 
the limit r — > 00, which corresponds to the limit u — > — 00, yield the same 
value, i.e., the total mass: 



lim Ex = lim £ 2 = M(-oo) = M . 



(51) 



20 



In fact it is proven in |H| that the total gravitational energies calculated 
out of triads related by a local SO (3) transformation, and that have the 
asymptotic behaviour given by (2), are the same. 

We consider next triads given by (25). The components of the torsion 
tensor are given by 



T (1) i2 = (--— +rd 1 (-)+^d 2 f)cos9cos ( j ) -(ld 1 f+^-f+-d 2 M)siri9cos(j ) , 
Vr r r 2r / \2 2r r J 

T(i)i3 — ( — I \- rd\(-)\sin9 sin(f) — —fcos9sin<j) , 



smc 



r r r J 2r 

^(1)23 = 2ccos9 sintp + (d 2 c — -f)sin9 sincfi , 

T (2)12 = (--—+rd 1 {-)+^d 2 f)ms9sin ( j ) -(ld 1 f+l^+-d 2 M)sin9 
Vrr r 2r / \2 2 r r J 

T( 2 )i3 = — (- H h rdi{-)\sin9 coscf) H f cos9cos<p , 

Vr r r / 2r 

^(2)23 = —2ccos9 coscf) — {d 2 c — -f)sin9 coscj) , 

T (3 )i2 = (—--- rd x {-) - ^-d 2 f)sin9 - (IdJ + ^-f+ -d 2 M]cos9 , 
Vrr r 2r / \2 2r r / 

2(3)13 = 1(3)23 = , (52) 

where <9i and d 2 denote partial derivatives with respect to r and 9, respec- 
tively. 

The calculation of (46) out of the components above does not pose any 
particular problem, except that the calculation is very long. Denoting by E 
the gravitational energy that follows from (25) and (52), we have 



E = E 1 - — 
2r 



7T 1 /"7T f)c f)c 

d9M 2 sin9 - - / cl9{— + 2ccotg9)—sin9 
4 Jo 09 '89 
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d9(— + Iccotgff) ccosO 



+ 7T 



8 7 



/ d9 (^E + 2c cotg8) 2 sin6{3cos 2 e - 1) 



(53) 



In the calculation above we have made use of the regularity condition c(u, 6) = 
for 9 = 0,7r. As expected, (53) also satisfies (51). 

The comparison of (49), (50) and (53) is crucial for the selection of triads. 
We are now in a position of obtaining numerical values for these expressions. 
With this purpose we make use of data based on the work of Saenz and 
Shapiro[|27j] f° r assigning values to Mq and u. In the latter reference it is 
discussed a model for stellar colapse of white-dwarves, with or without axial 
symmetry, rotating or nonrotating. In this model gravitational waves are 
produced during a burst of 10~ 3 s. Thus u may be taken to vary from 10~ 10 s 
to 10 _4 s. Based on the work of Saenz and Shapiro we consider a white- 
dwarf whose total mass is M = 1.4M , where M & is the solar mass. It is 
believed [^] that the maximum value for a white-dwarf is 1.4M . We will use 
geometrical unities (G=c=l) for the evaluation of the energies. The value of 
Mq in these unities is M & = 1.47664 x 10 5 cm. The distance r will be taken 
to vary from 10 18 cm to 10 13 cm. Finally, the parameters a and n assume the 
values 100 and —0.5, respectively, in proper unities. 

We have used the MAPLE V computer package to obtain Ei, E% and E. 
For the function F(fj,), which we now denote F(9), we take sin9, sin 2 '9, cos9 
and cos 2 9. The resulting values are listed in tables 1-4 for the particular 
value u = 10 _4 s. The numerical values were obtained with a precision of 20 
digits. 

We have verified that for all these functions, for all distances considered, 
E is always the minimum energy: 



Note that since Ex, E% and E are calculated at constant u = t — r, the radial 
dependence of these expression is given just by the ^ coefficient in (49), (50) 
and (53). Thus for constant u the intricacy of these expressions reside in 
the angular dependence. We have further verified numerically that the result 
above is also obtained for any value of u between 10 _10 s and 10 _4 s. The 
constants a and n must be chosen such that the absolute values of AM? 



E < E2 < E- 
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Table 1: F{0) = sinO 





r ~ 10 18 cm 


r ~ 10 15 cm 


r ~ 10 13 cm 


E l 


20653.912380952378 


20653.9123809523 


20653.9123809 


E2 


20653.912380952361 


20653.9123809350 


20653.9123792 


E 


20653.912380952344 


20653.9123809178 


20653.9123774 



Table 2: F(9) = sin 2 6 





r ~ 10 18 cm 


r ~ 10 15 cm 


r ~ 10 13 cm 


E 1 


20656.0287830687809 


20656.02878306878 


20656.0287830 


E 2 


20656.0287830687636 


20656.02878305149 


20656.0287813 


E 


20656.0287830687463 


20656.02878303421 


20656.0287696 



given by expressions (44) and (45) are not greater than M , otherwise the 
total gravitational energy will be negative. In any case, at present we do not 
know enough about axially-symmetric, nonrotating isolated sources in order 
to provide realistic values for these constants. 



VI. Discussion 

The result of the previous section constitutes a strong indication that 
the set of triads with asymptotic behaviour given by equation (2), and that 



Table 3: F(6) = cos9 





r ~ 10 18 cm 


r ~ 10 15 cm 


r ~ 10 13 cm 


E l 


20669.785396825396 


20669.785396825 


20669.78539683 


E 2 


20669.785396825379 


20669.785396808 


20669.78539509 


E 


20669.785396825362 


20669.785396790 


20669.78539336 
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Table 4: F(0) = cos 2 6 





r ~ 10 18 cm 


r ~ 10 15 cm 


r ~ 10 13 cm 


E 1 


20671.901798941799 


20671.9017989418 


20671.90179894 


El 


20671.901798941781 


20671.9017989245 


20671.90179721 


E 


20671.901798941764 


20671.9017989072 


20671.90179548 



satisfies the symmetry condition hij = hji, yields the minimum value for 
the gravitational energy that is computed from expression (9). In addition 
to the fact that the symmetrized triads are unique (namely, there does not 
exist a second set of triads that satisfy (2) and the symmetry condition), the 
present analysis indicates that the correct description of the gravitational 
field in terms of orthonormal triads, in the realm of the TEGR, is given by 
(24). The results described above amount to an interesting interplay between 
the energy properties of the space-time and its tetrad description. 

As long as we are interested in the dynamics of the gravitational field 
only, as described by the metric tensor, it is irrelevant which configuration 
of triads we adopt. However, if we consider the dynamics of spinor fields, 
such as the Dirac field, then the correct choice of triads is crucial. Recall 
that the theory defined by (5) and (6) was established under the asumption 
of the time gauge condition. The latter, together with (25), establishes the 
complete set of tetrad fields. And finally, if the detection of the emission of 
gravitational energy carried by gravitational waves is experimentally feasible, 
then the whole scheme developed here will play a relevant role. 

Bondi's radiating metric is valid only in the asymptotic region. There- 
fore we can do no better than determining the asymptotic behaviour of the 
triads. However in the more general case where the metric tensor is valid 
everywhere, except for singularities, we still have a prescription for assigning 
a unique set of triads to a given metric tensor restricted to the spacelike 
section. Such prescription is due to M0ller0, who called it "supplemen- 
tary conditions" . Although he established these conditions by still requiring 
de Donder relations for the metric tensor, we can dispense with the latter 
relations and ascribe generality to his proposal. M0ller suggested as supple- 
mentary conditions for the space-time tetrad field the weak field condition 
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e att « Vav + \h a » , (54) 

with h afl satisfying the symmetry condition h a ^ = h^ a . It differs from (2), 
which is restricted to spacelike sections only, in that the symmetry condition 
must be verified everywhere. Thus (54) is stronger than (2). Nevertheless we 
can require the set of triads to satisfy (54) in the general case. Every metric 
tensor for the spacelike section of a space-time admits a unique set of triads 
that satisfy a relation similar to (54), which is no longer a boundary condition 
and therefore it can be applied to space-times with arbitrary topology. We 
finally mention that the set of triads presented in ref. (jlfj, in the analysis 
of the irreducible mass of a rotating black hole, satisfies a relation similar to 
(54), but in the three-dimensional spacelike hypersurface. 
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